NUWC-NL  Technical  Report  10,237 
9  December  1992 


AD-A261  084 


Statistics  of  a  Whiteness  Measure 


Albert  H.  Nuttall 
Surface  ASW  Directorate 


Naval  Undersea  Warfare  Center  Detachment 

New  London,  Connecticut 

93-03628 

iililllli 

98  p  -  .. 


Approvad  for  public  r«leas«;  distribution  Is  unlimited. 


PREFACE 


This  research  was  conducted  under  NUWC  Project  Number  A70272 
Subproject  Number  RROONOO,  Selected  Statistical  Problems  in 
Acoustic  Signal  Processing,  Principal  Investigator  Dr.  Albert  H. 
Nuttall  (Code  302).  This  technical  report  was  prepared  with 
funds  provided  by  the  NUWC  In-House  Independent  Research  Program 
sponsored  by  the  Office  of  Naval  Research.  This  work  was  also 
sponsored  by  NUWC  Project  Number  A17653,  AN/BQR-22A  EC  15,  Task 
Assignment  Number  06U-93-7A432,  Sponsor  NAVSEA  06U23,  Project 
Manager  Evelyn  Hale. 


The  technical  reviewer  for  this  report  was  Alfredo  Edmonds 
(Code  2153). 


REVIEWED  AND  APPROVED:  9  DECEMBER  1992 


Donald  W.  Counsellor 
Director,  Surface  Antisulmiarine  Warfare 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No.  0704-0188 

PwMic  rcpoftin^  burden  for  this  coMeetiort  of  informAtion  is  estimeted  to  Avere^  i  hour  per  response,  including  the  time  for  reviewirtg  instructions,  searching  eiisting  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  renewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspea  of  this 
Election  of  information,  including  suggestions  for  r^ucing  this  burden  to  Washington  Headquarters  Services.  Oirectorate  for  information  Operations  and  Reports.  121S  iefferson 
Otvis  Highway.  Suite  1204.  Arlington.  VA  22202-4302.  and  to  the  Office  Qf  Management  and  Budget.  Paperwork  Reduction  Protect  (070441M).  Washington,  OC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blink)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  OATES  COVERED 

9  December  1992  Proqress 

4.  TITLE  AND  SUBTITLE 

Statistics  of  a  Whiteness  Measure 

S.  FUNDING  NUMBERS 

PE  61152N 

fi.  AUTHORCS) 

Albert  H.  Nuttall 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES} 

Naval  Undersea  Warfare  Center  Detachment 

New  London,  Connecticut  06320 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NUWe-NL  TR  10,237 

B.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AODRESS(ES) 

Chief  of  Naval  Research 

Office  of  Naval  Research 

Arlington,  VA  22217-5000 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

12a.  OISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release; 
distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Miximum  200  words) 

A  whiteness  measure  of  a  random  number  generator  is  defined 
as  the  sum  of  squares  of  all  the  off-zero  elements  of  the  sample 
covariance  function  of  a  finite  segment  of  data  of  length  K. 

The  mean  and  variance  of  this  whiteness  measure  are  evaluated 
exactly,  while  its  cumulative  and  exceedance  distribution 
functions  are  determined  by  simulations.  It  is  found  that  the 
variance  of  the  whiteness  measure  must  be  broken  into  the  two 
separate  cases  where  K  is  even  versus  K  is  odd.  This 
necessitates  the  analytic  evaluation  of  some  high-order  moments 
in  order  to  determine  the  variance  exactly. 

14.  SUBJEa  TERMS 

statistics  whiteness 

random  numbers  sample  covariance 

mean  variance 

IS.  NUMBER  OF  PAGES 

46 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 

IB.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

SAR 

NSN  7540-01  «280-S500 


Standard  ^orm  298  {Rev  2-89) 

P'McnbM  by  ANSI  Std 

2n-02 


TR  10237 


TABLE  OF  CONTENTS 

Page 


LIST  OF  ILLUSTRATIONS  ii 

LIST  OF  SYMBOLS  iii 

INTRODUCTION  1 

MEAN  AND  VARIANCE  OF  WHITENESS  MEASURE  3 

Mean  of  Whiteness  Measure  Wj^  4 

Variance  of  Whiteness  Measure  Wj^  5 

Special  Case  K  =  2  7 

Special  Case  K  =  3  7 

Special  Case  K  =  4  8 

General  Determination  of  Variance  of  Wj^  10 

PROBABILITY  DISTRIBUTIONS  OF  WHITENESS  MEASURE  15 

SUMMARY  25 

APPENDIX  A.  DERIVATION  OF  VARIANCE  OF  WHITENESS  MEASURE  Wj^  27 

APPENDIX  B.  PROGRAM  FOR  ESTIMATION  OF  DISTRIBUTIONS  OP  Wj^  37 

APPENDIX  C.  EVALUATION  OF  MOMENTS  DIRECTLY  FROM  MEASURED  39 

EXCEEDANCE  DISTRIBUTION 

REFERENCES  40 


DnC  QUAUTY  INSPFTTCD  3 


i 


TR  10237 


LIST  OF  ILLUSTRATIONS 

Figure  Page 


1. 

Variance  Vj^  for  Uniform  Random  Variables  {Xj^j 

13 

2. 

Variance  Vj^  for  Gaussian  Random  Variables  {Xj^l 

13 

3. 

Distributions 

of 

^2  for  Uniform  Random  Variables 

17 

4. 

Distributions 

of 

for  Uniform  Random  Variables 

17 

5. 

Distributions 

of 

W.  for  Uniform  Random  Variables 

4 

18 

6. 

Distributions 

of 

Wg  for  Uniform  Random  Variables 

18 

7. 

Distributions 

of 

for  Uniform  Random  Variables 

19 

8. 

Distributions 

of 

W22  for  Uniform  Random  Variables 

19 

9. 

Distributions 

of 

Wc.  for  Uniform  Random  Variables 

20 

10. 

Distributions 

of 

Wi28  for  Uniform  Random  Variables 

20 

11. 

Distributions 

of 

W2  for  Gaussian  Random  Variables 

21 

12. 

Distributions 

of 

for  Gaussian  Random  Variables 

21 

13. 

Distributions 

of 

W^  for  Gaussian  Random  Variables 

22 

14. 

Distributions 

of 

Wg  for  Gaussian  Random  Variables 

22 

15. 

Distributions 

of 

Wjg  for  Gaussian  Random  Variables 

23 

16. 

Distributions 

of 

W22  for  Gaussian  Random  Variables 

23 

17. 

Distributions 

of 

Wg^  for  Gaussian  Random  Variables 

24 

18. 

Distributions 

of 

Wi28  for  Gaussian  Random  Variables 

24 

Table  ' 

1.  Variance  of  Whiteness  Measure  Wj^  14 


11 


TR  10237 


LIST  OF  SYMBOLS 


K 

IID 
E(  ) 
F 


n 

Wr 

Var 

A#B,C,D 

A,B,C 

A,B,C 

P(x) 

D 

M 

u 

CDF(u) 

EDF(u) 

C 

E 


number  of  data  points,  (1) 
k-th  data  value,  (1) 
independent  identically  distributed 
expectation  of  random  variable,  (1) 
fourth  moment  of  random  variable  Xj^,  (1) 
n-th  delay,  (2) 

sample  covariance  at  delay  n,  (2) 
whiteness  measure,  (3) 
variance,  (10) 

variance  of  K-th  whiteness  measure  (10) 
unknown  constants  in  (10)  and  (35) 
sum  of  delayed  products  of  data,  (11) 
constants  for  K  odd,  (45) 
constants  for  K  even,  (48) 

probability  density  function  of  random  variable  x,  (51) 

constant  for  K  even  or  odd,  (55) 

size  of  fast  Fourier  transform,  (56) 

fast  Fourier  transform  of  data  (56) 

threshold  value,  (56) 

cumulative  distribution  function,  (57) 
exceedance  distribution  function,  (57) 
interval  probabilities,  (B-1) 
cumulative  distribution  function,  (B-3) 
exceedance  distribution  function,  (B-3) 
n-th  moment,  (B-8) 


iii/iv 

Reverse  Blank 


TR  10237 


STATISTICS  OF  A  WHITENESS  MEASURE 

INTRODUCTION 

When  a  random  number  generator  is  designed  to  yield  zero-mean 
independent  random  variables,  one  useful  test  of  its  validity  is 
afforded  by  its  sample  covariance  function.  This  quantity  would 
ideally  be  zero  for  all  delays  except  the  origin  value.  However, 
in  practice,  due  to  the  finite  length  of  data  generated  and  used 
to  test  the  generator,  the  sample  covariance  function  is  not 
identically  zero  but  fluctuates  about  zero.  A  measure  of  the 
whiteness  of  the  generator  is  afforded  by  the  sum  of  squares  of 
all  the  off-zero  elements  of  the  sample  covariance  function, 
relative  to  the  square  of  its  origin  value.  This  measure  was 
suggested  in  [1;  appendix  C]. 

In  this  report,  we  investigate  the  statistics  of  this 
whiteness  measure,  including  its  cumulative  and  exceedance 
distribution  functions  and  its  mean  and  variance.  Since  a 
sample  covariance  involves  products  of  data  values,  the  squared 
covariance  depends  on  fourth-order  products  of  the  data,  and  the 
variance  of  this  sample  quantity  involves  eighth-order  products 
of  the  data  under  various  delays.  It  is  this  latter  high-order 
product  which  greatly  complicates  the  statistical  analysis  and 
which  necessitates  a  roundabout  procedure  for  exact  evaluation  of 
the  variance  of  the  whiteness  measure.  The  probability 
distributions  of  this  measure  are  determined  by  simulation  for 
two  types  of  random  variables,  uniform  and  Gaussian. 
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MEAN  AND  VARIANCE  OF  WHITENESS  MEASURE 

Consider  real  data  sequence  Xq,  ...  of  K  data 

points  which  are  independent  and  identically  distributed  (IID) 
with  a  synunetric  probability  density  function  about  zero.  This 
zero-mean  sequence  will  have  all  odd-order  moments  equal  to  zero. 
Also,  assume  that  the  data  are  scaled  to  have  unit  variance  and  a 
fourth  moment  of  value  F;  that  is 

e(x^)  =  1  ,  e(xJ)  =  F  ,  for  0  i  k  <  K-1  ,  (1) 

where  E  denotes  the  expectation.  This  situation  includes  the 
uniform  random  number  generator  and  the  Gaussian  random  number 
generator,  for  example.  For  the  usual  uniform  random  variable 
distributed  over  (-^,is),  we  have  scaled  its  output  by  /12  for 
present  purposes  in  order  to  realize  variance  1.  Thus,  F  ®  1.8 
for  the  uniform  case,  while  F  =  3  for  Gaussian  numbers. 

The  sample  covariance  of  the  available  data  is  defined  as 

"n  '  R  ^  =‘k  =‘k-n  all  n  .  (2) 

Ideally,  we  might  like  to  have  sequence  {R^^l  equal  to  zero  for 
n  0,  However,  this  is  never  the  case,  although  the  (R^)  for 

n  0  are  much  smaller  than  Rq  when  K  is  large.  The  mean  value 

of  Rq  is  easily  seen  to  be  1,  by  reference  to  (1).  A  measure  of 
the  whiteness  of  data  sequence  Jxj^}  is  afforded  by  the  sum  of 
squares  of  all  the  off-zero  elements  of  sequence 

K— 1 

w„  H  FI  r2  =  2  ZZ  Rn  K  i  2  .  (3) 

^  W  ^  n=l  " 
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MEAN  OF  WHITENESS  MEASURE  Wj^ 

2 

The  mean  value  of  random  variable  R^^  follows  from  (2)  as 

-  ®(;;5  =‘k-n  *3  =‘j-J  - 

'  ri  ^  Cl  ^(xk  =‘k-n  =■:  *j-n)  ■ 

K  K  J 

Since  we  are  only  interested  in  values  of  n  >  0  according  to  ( 3 ) , 
the  expectation  in  (4)  is  nonzero  only  when  k  =  j;  here,  we  are 
utilizing  both  the  IID  and  the  zero-mean  properties  of 
Then,  (4)  becomes,  upon  use  of  (1), 

K—  1 

EfR^l  =  ^  HZ  1  =  —5^  for  1  <  n  1  K-1  .  (5) 

^  k=n 

(For  completeness,  E^RqI  =  (F  +  K  -  1)/K;  Variance(RQ)  =  (F-l)/K. 
Thus,  Rq  clusters  around  1  as  K  -♦  ®,  while  R^^  ->  0  as  K  -»  ®  for 
fixed  n  0.)  Use  of  result  (5)  in  (3)  yields  the  desired  mean 
value  of  whiteness  measure  Wj^  as 

K—  1 

^  HI  (K  -  n)  =  .  (6) 

K  n=l 

Notice  that  this  mean  value  is  independent  of  fourth-moment  F  and 
that  it  approaches.  1  as  K  ®.  Recall  that  E(Rq)  =  1  for 
comparison. 
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VARIANCE  OF  WHITENESS  MEASURE  Wj^ 

The  direct  evaluation  of  the  variance  of  random  variable 
in  (3)  would  require  a  very  tedious  procedure.  Whereas  the  mean 
evaluation  in  (4)  only  encountered  fourth-order  products  of 
delayed  versions  of  we  would  now  encounter  eighth-order 

products,  requiring  a  complicated  counting  procedure  to  account 
for  all  the  various  types  of  terms.  Specifically,  from  (2)  and 
(3),  we  have  whiteness  measure 


K-1  K-1  K-1 


Wk  =  -2  ^  L_  Xr  ^k-n  *j-n  ' 

K  n=l  k=n  j=n 


leading  to  mean  square  value 


.  2^  4  K-1  K-1  K-1  K-1  K-1  K-1 
'  K^  n=l  m=l  k=n  j=n  q=m  p=m 

E  x,  X,  X.  X.  X  x„  „  x„  x^  \ 

(  k  k-n  j  j-n  q  q-m  p  p-mj 


(7) 


(8) 


Not  only  would  this  eighth-order  average  have  to  be  evaluated  for 

all  possible  values  of  n,m,k,j,q,p,  but  the  sixth-order  summation 

would  then  have  to  be  conducted.  The  only  reasonable  case  that 

can  be  evaluated  from  (8)  is  that  for  the  term  proportional  to 
2 

F  .  It  is  obtained  only  for  the  special  choices  n  =  m  and 
k  =  j  =  q  =  p;  then  the  right-hand  side  of  (8)  reduces  to 


K-1 


(K  -  n) 


2(K  -  1)  „2 

.,3 


(9) 


Notice  that  moments  of  {Xj^j  above  the  fourth  need  not  be  known. 
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The  difficulty  of  attempting  to  evaluate  (8)  directly  forces 
us  to  attack  the  problem  from  a  different  aspect.  Specifically, 
we  adopt  a  shortcut  to  obtain,  exactly,  the  variance  of  whiteness 
measure  Wj^.  First,  observe  from  (8)  that  the  mean  square  value 

4 

of  Wj^  contains  a  denominator  of  K  .  Secondly,  it  has  been 
observed  from  simulations  that  the  variance  of  Wj^  goes  to  zero 
proportional  to  1/K  for  large  K.  Therefore,  the  form  of  the 
variance,  Vj^,  of  random  variable  must  be 


=  Var(Wj^) 


AK^+BK^+CK+D 


K 


(10) 


where  A,  B,  C,  D  are  unknown  constants.  In  order  to  determine 
these  four  constants,  we  will  evaluate,  exactly,  the  variance 
of  Wj^  for  a  sufficient  number  of  low-order  values  of  K,  and  then 
solve  the  four  simultaneous  linear  equations  yielded  by  (10). 

For  convenience,  we  define  the  sums 

♦n  =  |Z  ^k-n  ior  1  in<  K-1  .  (11) 

k=n 


Then 


Rn  =  ^  ♦n  1  <  n  i  K-1  , 


(12) 


as  seen  from  (2). 
form 


The  whiteness  measure  in  ( 3 )  then  takes  the 


0  K-1  „ 

Wk  =  ^  g 


for  K  i  2  . 


(13) 


For  K  =  1,  there  are  no  terms  in  the  sum,  yielding  =  0. 
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SPECIAL  CASE  K  =  2 

We  have,  from  (11)  and  (13), 


Therefore,  upon  use  of  the  IID  property  of  the  {Xj^l  and  (1), 

e(w2)  .  i  e(x<  xj)  .  {  .  (15) 

The  variance  of  W2  then  follows  as 

Vj  .  VarlWj)  -  e(w2]  -  =  i(F2  -  1)  ,  (16) 

where  we  used  ( 6 ) . 

SPECIAL  CASE  K  *  3 

The  procedure  for  the  remaining  cases  is  similar  to  that 
detailed  above  for  K  =  2;  therefore,  the  following  presentation 
will  be  abbreviated,  and  only  the  main  results  will  be  listed. 

We  have 


♦1  *  *1  *0 

+  X2  Xi  , 

^2  *2  *0  ' 

(17) 

W3 

2  f .2  .  ,21 

2f  2. 

.  »2  .  2  21 

(18) 

=  9lh  ^  ♦2I 

=  9lxi(Xo 

+  X2)  +  X2  XqJ  , 

.,2  4 

f  4.  ,  ,4 

,  4  4  , 

0  2 ,  ,2221 

(19) 

^^3  "  81 

1*11*0  *  *2' 

*2  ''O 

2  Xj(Xo  +  X2)  X2  XqI  . 

The  mean  value  of  (19)  is  given  by 
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^('^3)  “  +  6  +  F)  +  +  2(F  +  F)]  =  ^(3F^  +  IOf]  .  (20) 


Finally,  the  variance  of  is 


V  =  -If 

''3  81 1 


3F  +  lOF 


-  • 


(21) 


SPECIAL  CASE  K  =  4 

In  this  case,  we  have 

♦1  “  "^1  *0  ^2  ^1  *3  *2  '  ^2  “  *2  ^0  "^3  ^1  '  *3  '  *3  *0  ' 

(22) 

«4  “  uWl  *  *l  *  *1]  '  <“> 

64  wj  *  +  +2  ■*■  +3  2  ^2  ^  *3  ^  *2  *3  * 

The  mean  value  of  (24)  will  be  found  in  stages.  The  six 
components  of  (24)  have  the  following  average  values: 

^(♦3)  -  ^(='3  -S)  =  '  <25) 

^(♦3  ♦^l  *  "'O  ""3  ""l)^)  =  F  +  F  =  2F  ,  (26) 

e((|»3  =  E^Xj  Xq(x^  Xq  +  X2  Xj  +  X3  *2)^]  =  F+l+F  *  2F  +  1  , 

(27) 

e(+^)  =  e((X2  Xq  +  X3  Xj)^)  =  F^  +  6  +  F^  =  2F^  +  6  ,  (28) 
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®(*2  ♦?)  -  ^(<=‘2  =‘0  *  ==3  '='1  *0  +  *2  ^'l  *  ''S  ='2'^)  ' 

=  e([x2  *2  +  *2  +  2  [^2  ^2  +  ^2  ^2  ^ 

+  x2  x2  +  2  x2  x„  +  2  X3  Xj  Xj  x,,  +  2  X3  x^  Xj))  - 

=F+F+F+F+F+F+4=6F+4  ,  (29) 

♦1  “  *1<*0  *2^^  *3  *2  ■*■  ^  *3  *2  *1^*0  ^2^  ' 

44  444  222  2 

♦l  =  >‘l(*o  =‘2>  *  ’‘3  *2  ‘  =*3  ==2  =‘l<=‘0  ■"  *2> 

3  3  3  3 

+  4  X3  X2  Xj(Xq  +  X2)  +4  X3  X2  x^CXq  +  X2)  ,  (31) 

e(^j)  =  F(F  +  6  +  F)  +  F^  +  6(1  +  F)  »=  3F^  +  12F  +  6  .  (32) 

Combining  these  results  into  (24),  we  have  mean  square  value 

®K)  =  “) 

and  variance 

V4  =  3^(3f2  ^  ^ 

The  analytical  derivations  of  V^,  Vg,  Vy,  Vg  are  deferred 
to  appendix  A  due  to  their  lengthy  calculations  and  need  for  a 
shorthand  notation.  It  will  turn  out  that  we  also  need  all  of 
these  latter  results  when  we  find  the  constants  A,  B,  C,  D  in 
variance  expression  (10). 
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GENERAL  DETERMINATION  OF  VARIANCE  OF 

The  general  form  for  the  variance  Vj^  of  whiteness  measure 
is  given  by  (10)  for  arbitrary  K  and  is  repeated  below: 


=  Var(Wj^)  = 


AK^+BK^+CK+D 


(35) 


However,  analytic  determination  of  Vj^  for  K  =  2,  3,  4,  5,  6,  7,  8 
(see  appendix  A  also)  have  revealed  that  separate  forms  like  (35) 
must  be  employed  for  K  even  versus  K  odd.  That  is,  two  different 
sets  of  constants  A,  B,  C,  D  apply  in  the  even  versus  odd  cases 
of  K.  The  available  analytic  results  for  Vj^  (above  and  in 
appendix  A)  are  summarized  below: 


Vl  =  0 

(see 

i  the 

line  under  (13))  , 

(36) 

^^2  = 

-)  - 

(37) 

^3 

4 

~  81 

2 

3F‘^  + 

lOF  -  9]  , 

(38) 

^4 

1 

32 

3F^  + 

16F  -  7)  , 

(39) 

^5  = 

8  1 
■  6251 

5F^  + 

38F  -  11]  , 

(40) 

324  ^ 

L5F^  + 

144F  -  23]  , 

(41) 

4 

2401 

[21F^ 

+  246F  -  23]  , 

(42) 

^8 

1 

25( 

;(7r2 

+  96F  -  3]  . 

(43) 
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If  we  take  K  equal  to  the  odd  values  1,  3,  5,  7  in  (35)  and 
use  results  (36),  (38),  (40),  (42),  we  obtain  four  simultaneous 
linear  equations  for  the  constants  A,  B,  C,  D.  Their  solution 
leads  to  the  following  expression  for  the  variance  of  Wj^: 

+  B  K  +  c]  for  K  odd  ,  (44) 

K 

where 

A  *  4F  +  ^  ,  B  =  2F^  -4F--^,  C  =  -4F  +  8.  (45) 

When  (45)  is  substituted  into  (44),  the  variance  expression  can 
be  rearranged  in  terms  of  powers  of  F: 


V  =  K^, 1.1 

* 


KF^  +  2  K^-K-i]f  +  i(2K^-19K+12) ]  for  K  odd  .  (46) 


2 

The  F  term  here  confirms  (9),  as  anticipated. 

If  we  take  K  equal  to  the  even  values  2,  4,  6,  8  in  (35)  and 
use  results  (37),  (39),  (41),  (43),  we  obtain  four  different 
simultaneous  linear  equations  for  the  constants  A,  B,  C,  D. 

Their  solution  leads  to  the  following  expression  for  the  variance 
Vk  of  W^: 


U  =  +  B  K  +  c)  for  K  even  , 

K 


(47) 


where 


A  *  4F  +  ^  ,  B  =  2F^  -  8F  -  14  ,  C  = 


-  .  if 


(48) 
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When  (48)  is  substituted  into  (47),  the  variance  expression  can 
be  rearranged  in  terms  of  powers  of  F  according  to 


V  s 


(K-l)F^  +  2K(K-2)F  +  ■||2K^-2 lK+3l]  ]  for  K  even  .  (49; 


2 

Again,  the  F  dependence  in  (9)  is  confirmed. 

The  asymptotic  behavior  of  variance  Vj^  for  large  K  is  given 
by 


'K 


-  ("'■  *  3)  F  «  >' 


(50) 


for  both  K  odd  and  K  even.  This  is  due  to  the  fact  that  constant 

A  in  (35)  is  identical  for  the  odd  and  even  cases;  compare  (45) 

and  (48).  Thus,  whiteness  measure  W^^  tends  to  cluster  around  1 

as  K  -♦  *.  Recall  that  R^  1»  while  R„  0  for  fixed  n,  as  K  ®. 

u  n 

The  end  results  for  variance  Vj^  of  whiteness  measure  Wj^  are 
given  by  (44)  and  (47),  or  by  (46)  and  (49).  Plots  of  for  the 
uniform  random  variable  and  the  Gaussian  random  variable  {Xj^}  are 
displayed  in  figures  1  and  2,  respectively.  A  short  tabulation 
of  is  given  in  table  1  for  the  uniform,  Gaussian,  exponential, 
and  alternating  random  variables  Ixj^}.  The  probability  density 
functions  of  fxj^}  for  these  four  cases  are,  respectively. 


p^(x)  *  .5//3  for  |x|</3,  F*1.8;  (51) 

Pg(x)  =  (2n)"**  exp(-x^/2)  ,  F  *  3  ;  (52) 

Pe^*^  “  71  exp(-/2|x|)  ,  F  =  6  ;  (53) 

Pa(x)  =  ^  i(x-l)  +  7  6(x+l)  ,  F  =  1  .  (54) 
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0  4  %  11  1C  20  M  32. 

K 


Figure  1.  Variance  Vj^  for  Uniform  Random  Variables  {Xj^j 


0  -f  ?  /2  1C  zo  it  2?  32. 

K 

Figure  2.  Variance  Vj^  for  Gaussian  Random  Variables  jxj^j 
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A  short  table  of  the  variances  for  these  four  examples  is 
given  below.  For  the  alternating  example,  Xj^  =  ±1  and  F  =  1, 
whiteness  measure  W2  for  K  =  2  is  always  equal  to  1/2,  thereby 
leading  to  variance  V2  =  0.  The  smallest  possible  example  of  F 
is  1,  as  realized  in  the  alternating  random  variable  case. 


Table 

1 .  Variance 

of  Whiteness  Measure  W 

r 

K 

K 

Uniform 

Gaussian 

Exponential 

Alternating 

2 

.56 

2. 

8.75 

.0 

3 

.92444 

2.37037 

7.85185 

.19753 

4 

.985 

2.125 

6.15625 

.375 

5 

.94208 

1.8944 

5.0816 

.4096 

6 

.87901 

1.67901 

4.26235 

.41975 

7 

.81273 

1.50604 

3.68013 

.40650 

8 

.75188 

1.35938 

3.22266 

.39063 

16 

.45117 

.75586 

1.60986 

.25977 

32 

.24569 

.39722 

.79987 

.14771 

64 

.12804 

.20346 

.39808 

.07852 

128 

.06534 

.10295 

.19850 

.04045 

If  we  combine  (47)  with  the 

multiplied-out  version 

of  (44), 

the 

variance 

can  indeed  be  written  in  the  form  (35) 

for  all  K, 

where  the  constants  A,  B,  C  are  as  given  in  (48),  but  constant  D 


must  be  taken  according  to  the  two  different  values 


D 


0  for  K  evenj 
4 ( F  -  2 )  for  K  odd  J 


(55) 


Notice  that,  despite  (8)  involving  eighth-order  products,  nothing 
above  fourth-order  moment  F  of  is  required  in  these  results. 
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PROBABILITY  DISTRIBUTIONS  OF  WHITENESS  MEASURE 


The  direct  evaluation  of  whiteness  measure  according  to 
its  definition  (3)  in  conjunction  with  (2),  is  very  time 
consuming  for  large  K,  due  to  all  the  multiplications  required. 
An  attractive  alternative,  in  terms  of  fast  Fourier  transforms, 
was  derived  in  [1;  appendix  C]  and  is  employed  here;  the  program 
utilized  is  listed  in  appendix  B.  The  key  relation  relative  to 
(3)  is  [1;  (C-5)] 


W  S 


M-1 


M 


m=0 


IXml 


M-1 


m=0 


(56) 


where  M  is  the  size  of  the  fast  Fourier  transform  JX„J  of  data 
{Xj^l.  The  only  restriction  on  M  is  that  we  must  use  M  i  2K  -  1; 
then,  the  right-hand  side  of  (56)  is  independent  of  M.  (For 
K  *  1,  Xjjj  *=  Xq  for  0  i  m  i  M-1,  leading  to  Wj^  =  0,  as  noted  under 
(13).)  Again,  notice  that  the  whiteness  measure  depends  on 
fourth-order  products  of  the  data  or  its  transform. 

The  cumulative  distribution  function  (CDF)  and  exceedance 
distribution  function  (EDF)  of  whiteness  measure 

CDF(u)  *  Prob(Wj^  <  u)  ,  EDF(u)  *  Prob(Wj^  >  u)  ,  (57) 


for  the  case  where  data  fxj^)  is  uniformly  distributed  over  -/3,/3 
[see  (51)],  are  displayed  in  figures  3-10  for  K  =  2,  3,  4,  8, 
16,  32,  64,  128,  respectively.  These  results  were  determined  by 
using  at  least  one  million  trials  for  as  defined  in  (56).  The 
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exceedance  distribution  function  for  small  K  has  a  cusp  near  zero 
argument  which  disappears  for  larger  K.  However,  random  variable 
does  not  approach  Gaussian  as  K  increases;  rather,  as  shown  in 
figure  10  for  K  *  128,  the  right-hand  tail  appears  to  approach 
exponential  behavior.  For  a  bounded  random  variable,  |x^|  <  B, 
the  value  of  W„  is  bounded  according  to 


<  <■<  -  -  1)  _ 


(58) 


In  the  case  of  the  uniform  random  variable  Xj^,  where  B  *=  /3,  (58) 
yields  4.5  for  K  =  2,  10  for  K  =  3,  and  15.75  for  K  =  4. 

Although  the  mean  of  Wj^28  127/128  and  its  variance  is 

Vi28  *  .06534,  the  standard  deviation  of  Wj^28  0.256;  this 

leads  to  the  possibility  of  large  values  of  Wj28  on  occasion. 

For  example,  figure  10  shows  that  the  whiteness  measure  can  reach 
a  value  of  1.8  or  larger  about  1%  of  the  time.  If  a  candidate 
uniform  random  number  generator  has  probability  distributions  for 
Wj^  which  differ  significantly  from  figures  3  -  10,  it  is  suspect 
and  should  be  more  thoroughly  investigated  before  further  use. 

The  corresponding  cumulative  and  exceedance  distribution 
functions  of  the  whiteness  measure  for  a  Gaussian  random 
number  generator  (see  (52)]  are  displayed  in  figures  11-18  for 
K  *  2,  3,  4,  8,  16,  32,  64,  128,  respectively.  The  first 
observation  to  make  is  that  the  positive  tail  of  Wj^  can  now  reach 
much  larger  values  when  K  is  small.  However,  for  the  larger 
values  of  K,  the  probability  distributions  of  appear  to  be 
approaching  a  common  behavior,  regardless  of  the  distribution  of 
the  underlying  data  compare  figures  10  and  18  for  K  =  128. 


i.r  2  2r  3  3.  S’  4  4S- 

UL 


Distributions  of  W2  for  Uniform  Random  Variables 
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Figure  8.  Distributions  of  for  Uniform  Random  Variables 
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Figure  9.  Distributions  of  Wg^  for  Uniform  Random  Variables 
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Figure  10.  Distributions  of  for  Uniform  Random  Variables 
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SUMMARY 


The  statistics  of  a  whiteness  measure,  for  testing  a  random 
number  generator,  have  been  investigated  in  terms  of  the  mean, 
variance,  and  probability  distributions.  The  mean  and  variance 
results  are  exact  and  have  been  borne  out  by  numerous  simulations 
for  different  noise  sources  (X|^)  and  data  sizes  K.  These  results, 
for  whiteness  measure  defined  in  (3),  are  summarized  below: 


E,w^,  =  !4-i  ,  V,  .  var(«,,  .  A  ^  B  .  C  K  ^  D  ^ 


K 


where 


A  =  4F  +  I  ,  B  =  2F^  -  8F  -  14  ,  C  = 


-  2F^  +  ~  for  all  K  , 


while 


-I 


0  for  K  even 
4(F  -  2)  for  K  odd 


!■ 


(60) 


The  mean  of  whiteness  measure  Wj^  is  independent  of  fourth-order 
moment  F,  while  the  variance  of  depends  on  F,  but  not  on  sixth 
or  eighth-order  moments  of  data  That  is,  the  eighth-order 

product  encountered  in  the  general  mean-square  expression  (8) 
never  requires  knowledge  higher  than  fourth-order  for  its 
evaluation.  This  result  applies  for  a  symmetric  zero-mean 
probability  density  function  for  unit-variance  data 

The  cumulative  and  exceedance  probability  distributions  were 
determined  by  simulations  involving  more  than  one  million  trials 
each  and  therefore  have  good  reliability  approximately  down  to 
the  .0001  probability  level. 
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APPENDIX  A.  DERIVATION  OF  VARIANCE  OF  WHITENESS  MEASURE 

The  variances  Vj^  of  whiteness  measure  Wj^  for  K  *  2,  3,  4  were 
derived  in  (14)  -  (34)  in  the  main  text.  We  now  present  the 
derivations  for  the  remaining  cases,  K  =  5,  6,  1 ,  8,  that  are 
necessary  in  order  to  determine  V,.  for  all  K. 


SPECIAL  CASE  K  *  5 


♦l  '  *1  *0  *2  ^1  ^3  *2  *4  *3  '  *n  / 


4  ^0 


=  X2  Xq  +  X3  Xj  +  x^  X2  ,  ♦,  =  X,  Xo  +  X,  ,  (A-1) 


3  *0  *4  *1 


“5  "  ^  *3  *  +4]  ' 


(A-2) 


625  ,^2 
4  5 


J.  a.  ,  .2  .2  .  .2  .2  ,  ,  .2  .2  . 

♦1  ^  ^2  ^  *3  +4  2  ♦!  +2  *3  +  2  + 


2  2  2  2  2  2 
+  2  +  2  ^  2  +2 


(A-3) 


The  component  averages  required  are  developed  in  detail  as 
follows: 


E(tJ)  =  e(*J  xJ)  .  ,  (a-4) 

e(*4  fj)  ■=  e(xJ  x^(X3  Xp  +  X4  Xj)^]  .  F  +  F  =  2F  ,  (A-5) 

^(♦4  *2)  =  ®(*4  =‘^=‘3  ^'l  *  ='2'*0  *4''1  * 

-  e(x|  x2[x2  x2  +  x^lxo  +  +  2  X3  Xj  Xj(Xo  +  x^)]]  - 
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=  1  +  (F  +  F)  =  2F  +  1  , 


(A-6) 


=  e[xJ  XqCXjCXq  +  X2)  +  X3(X2  +  *4)]^)  “ 

^  *3(X2+x^)  +  2  X3  X3  (Xq+X2  )  (X2'*‘*4  )  I  ]  ” 


(F  +  1)  +  (1  +  F)  =  2F  +  2  , 


(A-7) 


£(♦3)  *  e([X3  Xq  +  X4  Xj]^)  «  F^  +  6  +  F^  «  2F^  +  6  ,  (A-8) 

e|^3  ^2!  “  *0  *4  1*3  *3  *2 ^ ^4  *0^^  )  * 

-  ^([*3  *0  *  *4  2  =‘4  ^3  ='1  =‘ol|='3  *1  *^='4  *  "'o'^ 

2  X3  X2  X^(X4  +  Xq)]]  =  F  +  (1+F)  +  F  +  (F+1)  =  4F  +  2  ,  (A-9) 

^(♦3  ^l)  “  ®(f*3  *0  *4  *3<*4  ' 


e([x2  x2  +  x2  x2  +  2  x^  X3  x^  x^]  [x3(X2  +  x^)  +  X3(X4  +  X2 ) 

+  2  X3  x^{X2  +  Xo^^*4  ■*■  *2^)  “ 

»  (1  +  F)  +  F(1  +  1)  +  F(1  +  1)  +  (F  +  1)  +  4  =  6F  +  6  ,  (A-10) 

£(♦2)  ”  ®(^*3  *1  *2^*4  ' 

«  F^  +  6(1  +  1)  +  F(F  +  6  +  F)  =  3F^+  6F  +  12,  (A-11) 
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e(^2  ^l]  *  e([X3  +  *2^*4  ■*■  *0  >  *3^*4  *2^^^) 

*  ^([*3  *1  "*"  *2^*4  "*”  *0^  "*”  ^  *3  ^2  ^1^*4  "*"  *0^]  ^ 

^  [^1^*2  "*”  *0^  ^  *3^*4  "*"  *2^  ^  ^  ^3  ^1^*2  ^  *0^^*4  "*"  *2^]]  * 

=  F(1  +  1)  +  F{1  +  1)  + 

+  e([x2  +  2  x^  Xq  +  x^][x^  +  2  x^  Xq  +  x^  Xq])  + 

+  e([xJ  +  2  x^  Xq  +  Xol[*4  *2  2  x^  x^  +  x^])  + 

+  4  e(x2(x^  +  Xq)(X2  +  Xq)(x^  +  X2))  = 

=4F+(F+1+F+F)+(F+F+1+F)+  4(1  +  1)  =  lOF  +  10, 

(A-12) 

^(♦l)  “  ®(f*l(*2  *0^  *  *3^*4  ■*■  *2^^^)  “ 

*  ^(*1^*2  *0^  "*"  ^  *1^*2  *0^  *3(^4  *2^  *3(^4  *2^  ]  * 

*  F{F  +  6  +  F)  +  6  e([x2  +  2  X2  Xq  +  ^q] [*4  +  2  x^  X2  +  x^])  + 

+  F(F  +  6  +  F)  *  4F^  +  12F  +  6(1  +  F  +  1  +  1)  « 

*  4F^  +  18F  +  18  .  (A-13) 
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Now,  we  combine  all  the  component  averages,  above,  to  obtain 
mean  square  value 

^(^s)  *  (A-14) 


and  variance 


V 


5 


— §— /  5F^ 

625'^' 


+  38F  -  11)  . 


{A-15) 


SPECIAL  CASE  K  *  6 

Now,  we  adopt  a  very  useful  shorthand  notation  to  handle  the 

rest  of  the  cases  of  interest.  For  example,  here,  ♦s  “  *5  *o 
2  2  2 

=  Xj  Xq,  which  IS  denoted  by  5500;  that  is,  the  superfluous  x 

2 

is  ignored  when  possible.  Also,  x^  X2  Xq  is  denoted  by  4220. 
With  this  background,  we  now  have 

♦5  *  5500  ,  «  4400+5511+2(5410)  , 

♦3  *  3300+4411+5522+2(4310+5320+5421)  , 

^2  *  2200+3311+4422+5533+2(3210+4220+5320+4321+5331+5432)  , 

-  1100+2211+3322+4433+5544+ 

+2(2110+3210+4310+5410+3221+4321+5421+4332+5432+5443)  .  (A-16) 

From  (13),  there  follows 

"6  -  ali  ♦?  -  TiK  +2  *  +5) 

n*i 
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and 


324  +  2{*lil  +  •••  +  ♦Jt^)  .  (A-18) 


We  also  abbreviate  the  following  ensemble  averages  as  follows 


E 


=  Tmn  . 


{A-19) 


Then,  there  follows,  in  a  straightforward  but  tedious  manner, 

T55  =  ,  T54  =  F+F  =  2F  ,  T53  »  F+l+F  =  2F+1  , 

T52  *  F+l+l+F  «  2F+2  ,  T51  =  F+l+l+l+F  =  2F+3  , 

T44  *  F^+6+F^  =  2F^+6  ,  T43  =  (F+F+1 )+( 1+F+F)  *  4F+2  , 

T42  *  (F+l+F+l)+{l+F+l+F)  =  4F+4  , 

T41  =  (F+l+l+F+F)+(F+F+l+l+F)+4  =  6F+8  , 

T33  »  3F^+4( 1+1+1 )+2{ 1+1+1)  »  3F^+18  , 

T32  -  (F+F+1+F)+(1+F+F+1)+(F+1+F+F)+4(1)  «  8F+8  , 

T31  «  (F+1+F+F+1)+(F+F+1+F+F)+(1+F+F+1+F)+4(1+1)  =  lOF+13  , 

T22  -  4F^+4(1+F+1+1+F+1)+2(1+F+1+1+F+1)  »  4F^+12F+24  , 

T21  =  (F+F+F+1+1 )+(F+F+F+F+l )+( l+F+F+F+F)+{ 1+1+F+F+F)+12  =  14F+18 
Til  =  5F^+(4+2) (F+l+l+l+F+l+l+F+l+F)  =  5F^+24F+36  .  {A-20) 

The  desired  average  is,  from  (A-18)  -  (A-20), 

324  E[wg)  «  15F^  +  144F  +  202  .  (A-21) 

The  variance  of  Wg  is  then 

''6  -  -  “)  • 


(A-22) 
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SPECIAL  CASE  K  =  7 

Continuing  in  the  fashion  established  above,  we  now  have 

=  6600  ,  =  5500+6611+2(6510)  , 

=  4400+5511+6622+2(5410+6420+6521)  , 

*  3300+4411+5522+6633+2(4310+5320+6330+5421+6431+6532)  , 

♦2  =  2200+3311+4422+5533+6644+ 
+2(3210+4220+5320+6420+4321+5331+6431+5432+6442+6543)  , 

=  1100+2211+3322+4433+5544+6655+2(2110+3210+4310+5410+ 
+6510+3221+4321+5421+6521+4332+5432+6532+5443+6543+6554)  .  (A-23) 

From  (13), 

«7  '  4f(*l  *  *2  "■•••*  ♦e) 

and  therefore 


MM  w2  -  *4  , 

4^  W7  -  ^1  + 


+  +  2[^l  +2  +5  ♦e)  •  (A-25) 


The  required  averages  are  as  follows: 

T66  =  ,  T65  =  F+F  *  2F  ,  T64  =  F+l+F  =  2F+1  , 

T63  =  F+l+l+F  =  2F+2  ,  T62  =  F+l+l+l+F  =  2F+3  , 

T61  =  F+l+l+l+l+F  =  2F+4  ,  T55  =  F^+F^+4+2  =  2F^+6  , 

T54  *  (F+F+1 )+( 1+F+F)  =  4F+2  ,  T53  *  (F+l+F+1 )+( 1+F+l+F)  =  4F+4  , 
T52  *  (F+l+l+F+l)+(l+F+l+l+F)  =  4F+6  , 

T51  =  (F+l+l+l+F+F)+(F+F+l+l+l+F)+4  =  6F+10  , 
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T44  =  3F^+4( 1+1+1 )+2( 1+1+1)  =  3F^+18  , 

T43  =  (F+F+l+l)+(l+F+F+l)+(l+l+F+F)  =  6F+6  , 

T42  *  (F+1+F+1+F)+(1+F+1+F+1)+(F+1+F+1+F)+4(1)  =  8F+11  , 

T41  *  (3F+3)+(4F+2)+(3F+3)+4( 1+1)  =  lOF+16  , 

T33  »  4F^+4(1+1+F+1+1+1)+2(1+1+F+1+1+1)  =  4F^+6F+30  , 

T32  =  (3F+2)+(3F+2)+(3F+2)+(3F+2)+4(l+l)  *  12F+16  , 

T31  «  (3F+3)+{4F+2)+(4F+2)+(3F+3)+4(l+l+l)  =  14F+22  , 

T22  *  5F^+(4+2) (l+F+l+l+l+F+l+l+F+l)  *  5F^+18F+42  , 

T21  *  2(3F+3)+3(4F+2)+4(l+l+l+l)  =  18F+28  , 

Til  *  6F^+(4+2) (F+l+l+l+l+F+l+l+l+F+l+l+F+l+F)  =  6F^+30F+60  . 

(A-26) 

The  average  of  interest  is,  from  {A-25)  and  (A-26), 

“  21F^  +  246F  +  418  ,  (A-27) 

leading  to  variance 

V7  =  24^11^^^^  *  (A-28) 
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SPECIAL  CASE  K  =  8 

This  is  the  last  case  that  we  need  to  evaluate.  We  now  have 

=  7700  ,  =  6600+7711+2(7610)  , 

=  5500+6611+7722+2(6510+7520+7621)  , 

=  4400+5511+6622+7733+2(5410+6420+7430+6521+7531+7632)  , 

^3  =  3300+4411+5522+6633+7744+ 
+2(4310+5320+6330+7430+5421+6431+7441+6532+7542+7643)  , 

+2  =  2200+3311+4422+5533+6644+7755+2(3210+4220+5320+6420+ 
+7520+4321+5331+6431+7531+5432+6442+7542+6543+7553+7654 )  , 

*  1100+2211+3322+4433+5544+6655+7766+ 
+2(2110+3210+4310+5410+6510+7610+3221+4321+5421+6521+7621+ 
+4332+5432+6532+7632+5443+6543+7643+6554+7654+7665)  .  (A-29) 

From  (13)  again, 

«8  '  6?(+?  ^  +2  *  +7)  ' 

giving 

1024  Wg  =  +  •••  +  +  2[tj  +2  "^  •••  +  ♦g  ♦7)  .  (A-31) 

The  averages  needed  are  listed  below. 

T77  =  ,  T76  =  2F  ,  T75  =  2F+1  ,  T74  =  2F+2  , 

T73  =  2F+3  ,  T72  *  2F+4  ,  T71  *  2F+5  , 
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T66  *=  F^+F^+4+2  =  2F^+6  ,  T65  =  ( F+F+1 )  +  ( 1+F+F )  =  4F+2  , 

T64  =  (F+l+F+l)+(l+F+l+F)  =  4F+4  , 

T63  =  (F+l+li-F+l)  +  (l+F+l+l+F)  =  4F+6  , 

T62  =  (F+l+l+l+F+l)+(l+F+l+l+l+F)  =  4F+8  , 

T61  =  (3F+4)+(3F+4)+4  =  6F+12  , 

T55  =  3F^+4( l+l+l)+2{ 1+1+1)  =  3F^+18  , 

T54  *  (F+F+l+l)+( l+F+F+l)+( 1+1+F+F)  =  6F+6  , 

T53  =  3(2F+3)  =  6F+9  , 

T52  =  {3F+3)+(2F+4)+(3F+3)+4{ 1)  =  8F+14  , 

T51  =  (3F+4)+(4F+3)+(3F+4)+4(l+l)  =  lOF+19  , 

T44  =  4F^+4{6)+2(6)  =  4F^+36  , 

T43  =  (3F+2)+{2F+3)+(2F+3)+(3F+2)+4(l)  =  lOF+14  , 

T42  =  4(3F+3)+4(l+l)  =  12F+20  , 

T41  =  2(3F+4)+2(4F+3)+4(l+l+l)  =  14F+26  , 

T33  =  5F^+4(2F+8)+2(2F+8)  =  5F^+12F+48  , 

T32  =  4(3F+3)+(4F+2)+4(l+l+l)  =  16F+26  , 

T31  =  2(3F+4)+3(4F+3)+4(l+l+l+l)  =  18F+33  , 

T22  =  6F^+4(4F+11)+2(4F+11)  =  6F^+24F+66  , 

T21  =  2(3F+4)+4(4F+3)+4(l+l+l+l+l)  =  22F+40  , 

Til  =  7F^+4(6F+15)+2(6F+15)  =  7F^+36F+90  .  (A-32) 

The  desired  average  is  therefore 

1024  e[Wq)  =  28F^  +  384F  +  772  ,  (A-33) 

giving  variance 

''8  ' 
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APPENDIX  B.  PROGRAM  FOR  ESTIMATION  OF  DISTRIBUTIONS  OF  W 


10 

T=1E6 

!  NUMBER  OF  TRIALS  "NUWC  TR10237" 

20 

K=32 

!  NUMBER  OF  DATA  POINTS,  ARBITRARY 

30 

M=64 

1  FFT  SIZE,  M  >s  2K-1,  POWER 

OF  2 

40 

L=1 1000 

!  NUMBER  OF  LEVELS  FOR  DISTRIBUTION 

50 

Dus , 001 

!  INCREMENT  IN  U 

60 

Grsl000 

!  GRID  SPACING 

70 

PRINTER  IS  PRT 

80 

PRINT  "K  =";K;" 

T  =";T}“  Du  =";Du;"  UNIFORM" 

90 

PRINTER  IS  CRT 

100 

DOUBLE  T,K,M,L,M1 

,Ms,M2,Kl,Ts,Ks  !  INTEGERS,  NOT 

DP 

110 

DIM  CosC512>,X(2048) , Y<2048), V<30000) 

120 

MlsM-1 

130 

REDIM  Cos<0:M/4), 

X<0;M1),Y<0:M1>,V<0!L) 

140 

fl=2.*PI/M 

150 

FOR  Ms=0  TO  M/4 

160 

Cos<Ms)=COS<fl*Ms) 

!  QUARTER-COSINE  TABLE  IN  Cos 

(*) 

170 

NEXT  Ms 

180 

M2*M/2 

190 

K1=K-1 

200 

Tl*! . /T 

210 

Fsl2. /<K*N) 

!  UNIT-VARIANCE  UNIFORM 

220 

F=F4F 

!  RANDOM  VARIABLES  {x<subk>> 

230 

MusKl/K 

!  EXACT  MEAN  OF  WK 

240 

Mul=Var»0. 

250 

TasTIMEDflTE 

260 

FOR  Ts=l  TO  T 

270 

FOR  Ks=0  TO  K1 

280 

XCKs)=RND-. 5 

!  ZERO  MEAN 

290 

Y<Ks)=0. 

!  REAL  INPUT 

300 

NEXT  Ks 

310 

FOR  Ks=K  TO  Ml 

320 

X<Ks)=Y<Ks)=0. 

330 

NEXT  Ks 

340 

CALL  Ffll4<M,Cos':*),X<*),Y<*)) 

350 

S2=S4=0. 

360 

FOR  Ms=l  TO  M2-1 

!  ZERO  TO  NYQUIST 

370 

X=X<Ms) 

380 

Y=Y(Ms> 

390 

fi=X*X+Y*Y 

400 

S2=S2+fl 

410 

S4=S4+fl*fl 

420 

NEXT  Ms 

430 

X»X<0) 

440 

fl«X(M2) 

450 

XsX*X 

460 

470 

S2=X+fl+2. *S2 

480 

S4=X*X+fl»fl+2. *S4 

490 

W=F*<M*S4-S2*S2) 

!  WHITENESS  MEASURE  WK 
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see  Mui-Muit-u 

Sie  Var=Var+(W-Mu>*<W-Nu) !  USE  KNOWN  MEfiN  Mu 
52e  MssINKW/Dw) 

53e  MsaMlN<Ms,L> 

54e  V<M£)*V<Ms)+T1  !  INCREMENTAL  PROBABILITIES 

SSe  NEXT  Ts 

see  Tb=TIMEDRTE 

S7e  PRINTER  IS  PRT 

see  PRINT  <Tb-Ta).''3be0;  "  HOURS" 

see  PRINT 

see  Mul=Mul/'T  !  ESTIMATED  MEAN  OF  WK 

6ie  Var=Var/T  •  ESTIMATED  VARIANCE  OF  WK 

62e  PRINT  "Mul  =";Mul;“  Mu  =“;Mu 

63e  PRINT  "Var  =";Var 

64e  PRINT 

6Se  PLOTTER  IS  "GRAPHICS" 

66e  GRAPHICS  ON 

67Q  WINDOW  e,L,-6,e 

see  LINE  TYPE  3 

ese  GRID  Gr, 1 

7ee  LINE  TYPE  1 

7ie  c=e. 

72e  FOR  Ms=e  TO  L-1 

730  C=C+V<M£)  !  CDF  OF  WHITENESS  MEASURE  WK 

740  IF  C>0.  THEN  760 

7S0  GOTO  770 

760  PLOT  Ms+1,LGT<C) 

770  NEXT  Ms 

780  PENUP 

790  E«S1«S2*0. 

800  FOR  Ms=L  TO  1  STEP  -1 

810  E»E+V<Ms)  !  EDF  OF  WHITENESS  MEASURE  WK 

820  S1«S1-^E 

830  S2=S2+Sl 

840  IF  E>0.  THEN  860 

8S0  GOTO  870 

860  PLOT  Ms.LGTCE) 

870  NEXT  Ms 

880  PLOT  0,0 

890  PENUP 

900  Mul=Dw*< . S+Sl )  !  ESTIMATED  MEAN  OF  WK 

910  Mu2=2. *Dw*Du*S2  !  SEE  APPENDIX  C 

920  PRINT  "Mul  =";Mul;"  Mu  =";Mu 

930  PRINT  "Var  a";Mu2-Mu*Mu  !  ESTIMATED  VARIANCE  OF  WK 

940  PRINT 

9S0  PRINTER  IS  CRT 

960  PAUSE 

970  END 

980  ! 

990  SUB  Fftl4<D0UELE  N,REAL  Cos < » > , X< # > , Y< * ) )  !  N<=2  14=16384;  0  SUBS 
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APPENDIX  C.  EVALUATION  OF  MOMENTS  DIRECTLY 
FROM  MEASURED  EXCEEDANCE  DISTRIBUTION 

Let  X  be  a  positive  random  variable  with  probability  density 
function  p,  cumulative  distribution  function  (CDF)  C,  and 
exceedance  distribution  function  (EDF)  E.  Let  the  measurements 
of  these  distributions  be  the  interval  probabilities 


E(A)  *  1-Vq  ,  E(2A)  =  I-Vq-V^  ,  E(3A)  =  1-Vq-V^-V2 


(C-7) 
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1 

The  first  two  moments  of  random  variable  x  can  be  developed 


as 


J  dx  X  p(x)  =  J  dx  E(x) 
0  0 


(C-8) 


and 


^2  *  J  dx  x^  P(x)  ®  2  J  dx  X  E(x)  = 
0  0 


2  A  E(nA)  .  (C-9) 

n=l 


These  results  can  be  rapidly  evaluated  by  recursion.  For 
E( (N+1)A)  =  0,  use 


E=S1=S2=0. 

FOR  Ns=N  TO  1  STEP  -1 
E=E+V(Ns) 

S1=S1+E 
S2=S2+S1 
NEXT  Ns 

Mul=Delta*( .5+Sl) 

Mu2=2.*Delta*Delta*S2  (C-IO) 
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